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We showed that quantum mechanical superposition can sustain in the process of ion transfer in 
protein membrane for a substantial period in spite of the presence of the interactions with environ- 
mental modes of molecular vibration. The spectral temperature, as defined in quantum thermody- 
namical framework plays a significant role in maintaining the coherence. The ratio of decoherence 
time and dwell time has been calculated, which can be directly related to the degree of coherence. 
The results shead new light to build quantum information system of entangled ionic states in the 
voltage gated biological channels. 
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PACS numbers: 03.65.Xp, 03.65.Yz, 87.10. Ca 



The approaches towards dynamics of ion transport in 
protein membranes (ion channels) are generally consid- 
ered to be classical and based on molecular dynamics or 
Brownian dynamics. Recently, one of the present au- 
thors showed that IJ the dynamics of K-ion channel can 
be explained using nonlinear Schrodinger equation which 
is compatible with the results of MacKinnon's experi- 
mental observation 0. Here, the two K-ions may form 
an entangled state within selectivity filter during a finite 
period of time. The temperature within the channel is 
generally considered to be high enough to destroy the 
coherence within very short period. Moreover, molecu- 
lar modes of the protein environment induce dynamical 
decoherence, which destroys the quantum mechanical su- 
perposition of states in a very short period of time. So 
quantum mechanical approach in these area of research 
was mostly speculative and far from experimental realiza- 
tion. But in recent times, some experimental demonstra- 
tion of the presence of quantum coherence in the pro- 
cess of photosynthetic energy transfer [1, 01 lead us to 
reconsider the theoretical approach and understanding. 
Here we are concerned with the question that whether 
and under what condition a sustainable quantum super- 
position is achievable in the process of transfer of ions 
through biological channels. Now it is quite a practical 
argument that the quantum state of the traversing ion 



is strongly coupled with the molecular vibrational modes 
of the protein environment and hence fast decoherence 
is almost absolutely unavoidable. But it is to be noted 
that the traversal time of the ion through the membrane 
is also quite small and it is the ratio of the time scale 
of decoherence with this traversal time that plays an im- 
portant role in understanding the maintenance of coher- 
ence. If the decoherence time is larger than the traversal 
time of the ion, then the quantum superposition of the 
ionic states is sustainable enough for the traversing en- 
tity within the period of ionic transfer. Here we also 
need to consider the effect of temperature in estimating 
the decoherence time. With the recent developments of 
quantum thermodynamics , a new concept of tempera- 
ture (known as spectral temperature) for micro-states at 
non-equilibrium condition has been proposed, instead of 
the usual concept of thermodynamic temperature. Since 
ion transport through protein membrane is essentially a 
non-equilibrium phenomena, we propose that the spec- 
tral temperature plays important role in understanding 
the coherence in the channel dynamics. 

Considering the master equation of the density opera- 
tor for a certain quantum system, the decoherence time 
[dl can be written as 
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where 7 is the relaxation (dissipation) parameter, T is 
the thermodynamic temperature and Aa; — x — x' 1% the 
spatial shift of the particle. Now there are numerous def- 
inition of traversal time, among which phase time and 
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dwell time are generally accepted by the community 0- 
Phase time is equated to dwell time with an additive 
self-interference term and argued to be the time interval 
between the energy storage and release in the barrier re- 
gion [S- So for non-dispersive barrier, where the self- 
interference term vanishes, the dwell time can be readily 
interpreted as the life time of energy stora ge a nd release 
in the barrier region. In a previous work 12|, we have 
calculated the weak value of dwell time as 



icothf^-^l 



7 



A 2 / 



(0.2) 



where tm is the measurement time. 
The protein membrane can be assumed as an array of 
molecules with certain modes of vibration. We have men- 
tioned earlier that the quantum states of the traversing 
entity is interactively coupled to the molecular vibration 
modes of the protein membrane. The traversing quan- 
tum entity loses energy to the vibrational modes, due 
to the presence of this environmental coupling. So the 
protein membrane can be interpreted as sort of capac- 
itive system, which can store and release energy with 
the dwell time as it's lifetime of energy storage. In this 
case the time interval between the opening and closing 
of the channel gates can be taken as the measurement 
time (tm ) ■ N ow we consider a double well potential of 
the form 13 1 
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where A and B are dimensionless constants. For the par- 
ticular case of asymmetric double well, A — 14 and 
B = 45. The bistable nature of the potential is use- 
ful in practical situation, since the ionic transfer in the 
channel can be interpreted as tunneling between the two 
stability regions situated at xq and X2 (see fig.l). 
It is also very important to reconsider the concept of 
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FIG. 1: V(x) vs. X with parameters A = 14 and B — 45 

temperature in this aspect. Temperature, as we know 
from the context of thermodynamics, is a property of 
equilibrium. The ensemble mean of the kinetic energy 



equals to Boltzman constant times the temperature. So 
assuming ergodicity, temperature is defined as the time 
averaged kinetic energy. But the process of ion trans- 
fer through channels is a dynamical process interactive 
to the protein environment and subject to energy ex- 
change with the environment. This is certainly not a 
situation for thermodynamic equilibrium and hence the 
usual concept of thermodynamic temperature may not 
be suitable. Here we introduce the concept of spectral 
temperature, originally formulated by Gemmer et.al jlj. 
It is defined as a function of the microstates, to include 
the non-equilibrium properties of the system. It is for- 
mulated as a function of energy occupation probability 
of the different states of a certain quantum system. This 
temperature evolves with the evolution of the probability 
of occurrence. The inverse of the spectral temperature is 
defined as Q 
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where Pj is the probability of finding the particle within 
an energy compartment with mean energy Ei , having the 
degree of degeneracy 4>i. This definition depends on the 
energy probability distribution and of course, the spec- 
trum of the concerning system. So it cannot change in 
time for an isolated system and defined independent of 
the fact that whether the system is in equilibrium or not. 
In this definition, the association of quantum probability, 
which evolves in time, gives temperature a evolving fea- 
ture representing the dynamical situation between two 
successive equilibrium. Now here we are approximating 
the bistable potential as a two-state system having only 
the ground states of the asymmetric wells. So for our 
case of non-degenerate two-state system, the expression 
of the spectral temperature reduces to 
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Where Pq and Pi are the probabilities corresponding to 
the ground states of the lower and higher well respec- 
tively. Since Pq is the probability corresponding to the 
lowest state, it should not decay. Now using the decay 
probability of the higher state within the time interval 
equal to the dwell time, we find 
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Now if we consider the energy loss by the particle travers- 
ing from the higher barrier to the lower one in terms of 
the usual kinetic temperature (T^) as Ei — Eq = ^KBTk, 
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we get 
T 
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FIG. 2: Tspec/Tk vs 7tm. Here we keep 7 as a constant and 
basically study the variation with increasing tm- 

From eqn (10. 7p . under the condition tm ^ we find 
that Tspec — Tk- ie. the spectral temperature is almost 
equal to the usual kinetic temperature, if the time inter- 
val of gate opening and closing is very greater than the 
dissipation time scale. If the gate opening and closing 
mechanism is slow enough to be considered as a quasi- 
static process, the spectral temperature remains very 
close to it's kinetic counterpart. Now it should also be 
noted that though the process is sort of quasi-static one, 
it should not be considered as a reversible process, be- 
cause the coupling to the molecular vibrational modes of 
the protein environment ensures some generation of dis- 
sipative entropy. 

In case of a double well potential given by (|0.3|) . using 
(I03]),(I021) and dOJl) we get 
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where eo = -Bi — Eq is the asymmetry energy of the po- 
tential and w — is the separation length between the 
wells. For the quasi-static condition tm ^ - 
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So we find that in the quasi-static region, whether the de- 
coherence time scale is larger than the dwell time depends 
on the mass of the traversing particle, the length scale of 
the channel and the asymmetry energy which is basi- 
cally the energy lost by the particle during the process 
of traversal. The ratio of the time scales is inversely pro- 
portional to all the mentioned parameters, i.e. greater 
the inertia of the traversing particle, stronger the inter- 
action (hence greater energy loss) and larger the traver- 
sal length, faster shall be the decoherence and hence the 
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FIG. 3: Tdec/TD VS 'jTM ■ Here also we keep 7 as a constant 
and basically study the variation with increasing tm- Here 

F = — , —2—. As with increment of tm, the process becomes 

quasi-static, we see that the ratio of the two timescales also 
reach a stable value 



situation will be more and more classical. But if the pa- 
rameter can be chosen in such a way in which the process 
of decoherence is slower than the process of ionic transfer 
(i.e. Tdec/TD > 1), then quantum superposition will be 
sustainable enough within the time period of ionic trans- 
fer. 

Now we turn our attention to the "degree of coher- 
ence" for such cases of ionic transfer through protein 
membranes and try to establish a relation of the afore 
mentioned ratio of decoherence-dwell time scales with 
it. Convenient way to discuss coherence in quantitative 
terms can be done through the introduction of normal- 
ized form of correlation functions such as "degree of co- 
herence" 
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where ^ = x/a. For our case of double well system ap- 
proximated as two-state system, (|0.10l) is reduced to 
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For completely coherent situation, the degree of coher- 
ence is 1 and for completely decohered case it is 0. Gen- 
erally the value should lie in between. If we approximate 
the potential around the left well at ^ = as a harmonic 
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potential, we find that the wavefunction can be estimated 
as [H 
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where u 
as 



Expressing the asymmetry energy 
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we find that for a double well potential {A — 14 and 
B = 45) 
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So for the interval of dwell time given by equation 
we find that the degree of coherence 
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For the quasi-static condition tm ^1/7 
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From equation (|0.9p and (|0.18p . we can also establish a 
relation 
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From this relation we see that for completely coherent 
and decoherent situation, the decoherence- dwell time ra- 
tio is infinity and zero respectively and in general situa- 
tion it lies in between. So this time-scale ratio gives us a 
certain measure of coherence. As the value of this ratio 
gets bigger, the "quantumness" of the traversing entity 
increases. 
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FIG. 4: Tdec/ro vs. g 



Depending on the above analysis, we suggest that the ion 



selectivity filter may exhibit quantum coherence which 
can play crucial role in the process of selectivity and con- 
duction of specific ions in biological membranes. For the 
time scales shorter than that of decoherence time, quan- 
tum coherence can be expected to sustain and have vital 
importance in the dynamics, despite the presence of in- 
teractive protein environment. Our analysis shows that 
for a sort of quasi-static situation, where the gate open- 
ing and closing mechanism is slower than the relaxation 
(dissipation) time scale, the decoherence-dwell time ra- 
tio reaches a static value and can also be greater than 
unity depending on the mass, energy and length param- 
eters. In such situations, coherent phenomena like en- 
tanglement can be of vital importance in understanding 
the mechanism of selectivity and transport. In case of 
K+ filter, there exists two energetically almost degener- 
ate binding states, commonly referred as (1,3) and (2, 4) 
states [l5-17|. Presence of quantum superposition may 
lead us to explain the transport phenomena in terms of 
quantum mechanical tunneling between these two states. 
The interplay between quantum coherence and environ- 
mental noise induced dephasing may also be of funda- 
mental importance. Progress in the atomic spectroscopy 
of membrane proteins indicate that protein membrane 
organization may carry a certain coding potency, im- 
plying quantum entanglement within ion channels 
Il8j . Increasing number of studies are indicating towards 
the probabilistic nature of ion channel gating mechanism 
0, [l§|. In the light of these researches, we conclude that 
it may be relevant to build certain model of quantum 
information system driven by the entangled ionic states 
in the voltage gated selectivity filters, which can provide 
necessary inferences into the biological ion channel dy- 
namics. 
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